APPENDIX 1 



Program for Baselining. Normalizing. Interpolating Then 
Calculating Spectral Overlap Integrals 
5 C This program has a non-standard DO WHILE loop 

INTEGER NPTS, NMAX, ROWS, ITER 
INTEGER EOF1, FLERR1, FLERR2 
INTEGER EOF2, FLERR3, FLERR4 
1 0 INTEGER EOF3, FLERR5, FLERR6 

INTEGER FLERR7 
INTEGER ij 

CHARACTER*30 fhamel, fhame2, fhame3, fhame4 
CHARACTER*30 fhame5, fhame6, fhame7 
N 15 P ARAMETER(NM AX=3 500, LAMDA=601) 

2 REAL x,xx 1 (NMAX) ,yyl (NMAX),INTERV 1 

H REAL xx2 (NMAX) ,yy2 (NMAX), INTERV2 

Jj REAL xx3 (NMAX), yy3 (NMAX),INTERV3 

Cf 20 REAL yil (NMAX), yi2 (NMAX), yi3 (NMAX), yc (NMAX) 

O REAL area 

■ CHARACTER* 1 SUBSTR, INITAR, LIGHT, INTMED 

^ FLERR1=0 
m25 FLERR2=0 
5 FLERR3=0 
H FLERR4=0 

FLERR5=0 

FLERR6=0 
30 FLERR7=0 

EOF 1=0 

EOF2=0 

EOF3=0 

INTERV1=0 
35 INTERV2=0 

INTERV3=0 

area=0 

write(*,*) 'Do you wish to output intermediate files? (Y/N)' 
40 read(*'(A)') INTMED 

write(*,*) 'Do you wish to process a substrate file? (Y/N) 5 
read(*,'(A)') SUBSTR 

45 IF ((SUBSTR.EQ. *Y') .OR. (SUBSTREQ. 'y ')) THEN 
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ITER=0 

do5ITER=l,NMAX 

xxl (ITER)=0 

yyl (ITER)=0 

yii (ITER)=0 
5 continue 

write(*,*) 'Enter the name of the input substrate file:' 
read(*,' (A)') fhamel 

open 

(TJNIT=ll,FILE=mamel,STATUS='OLD',IOSTAT=FLERRl,E 
RR=101) 

ROWS=0 

do while ((EOF1.EQ.0) .AND. (ROWS.LT.NMAX)) 
ROWS=ROWS+l 

Read (1 l,*,IOSTAT=EOFl) xxl (ROWS), yyl (ROWS) 

end do 

close (UNIT=11) 
NPTS=0 

IF (EOF1.NE.0) THEN 
NPTS=ROWS-l 

write(*,'(I4,lX,A12)') NPTS, 'points read.' 

ELSE 

NPTS=ROWS 

write(*,'(A28, 14, A12)') 'Too many data points! First', 
+ NMAX,' points read. . . ' 

END IF 

call baseln(yyl,NPTS) 

IF ((INTMED.EQ. ' Y') .OR. (INTMED.EQ. 'y')) THEN 

write(*,*) 'Enter the name of the output substrate file:' 

read(*,' (A)') fhame2 

open 

(UNIT=12,FILE=fname2,STATUS='NEW',IOSTAT=FLERR2, 
ERR=102) 

write(*,*) 'Writing data. . . ' 

END IF 
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do 12 i=l, LAMDA 

x= (i-l)+200 

call locate(xxl,NPTS,xj) 

if (O'.eq.O). OR. (j.eq.NPTS)) then 
INTERV1=0 

else 

INTERV1= ((yyl(j+l)-yyl(j)) / (xxl(j+l) - xxl(j))) * 
(x-xxlO)) 

+ +yylO) 
end if 

yil (i)=INTERVl 

IF ((INTMED.EQ. ' Y') .OR. (INTMED.EQ. 'y ')) THEN 

if (O'.eq.O) .OR. (j.eq.NPTS)) then 

GO TO 12 

else 

write(12,'(lx, f7.2,i6,3fl2.2)')xj,xxl(j),xxl0+l), 
+ INTERV1 

endif 

END IF 

12 continue 

IF ((INTMED.EQ. ' Y') .OR. (INTMED.EQ. 'y')) THEN 
close (UNIT=12) 

ENDIF 

ELSE 

ITER=0 

do 14 ITER=1, LAMDA 
yil (ITER)=1 
14 continue 

ENDIF 
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write(*,*) 'Do you wish to process an initiator file? (Y/N)' 
read(*,'(A)') INITAR 

IF ((INITAR.EQ.'Y') .OR. (INITAR.EQ.V)) THEN 
ITER=0 

dol5ITER=l,NMAX 

xx2 (ITER) =0 

yy2 (ITER)=0 

yi2 (ITER) =0 
1 5 continue 

write(*,*) 'Enter the name of the initiator file:' 
read(*,'(A)') mame3 

open 

0LJNIT=13,FILE=marne3 ) STATUS='OLD',IOSTAT=FLERR3 J 
ERR=103) 

ROWS=0 

do while ((EOF2.EQ.0) .AND. (ROWS.LT.NMAX)) 
ROWS=ROWS+l 

read(13,*,IOSTAT=EOF2) xx2 (ROWS), yy2 (ROWS) 

end do 

close(UNIT=13) 
NPTS=0 

IF (EOF2.NE.0) THEN 
NPTS=ROWS-l 

write(*,l(14,lX,A12)') NPTS,'points read.' 

ELSE 

NPTS=ROWS 

write (*,'(A28,I4,A12)') 'Too many data points! First 
+ NMAX, ' points read 

END IF 

call baseln (yy2,NPTS) 

IF ((INTMED.EQ.'Y').OR. (INTMED.EQ.'y')) THEN 

write(*,*) 'Enter the name of the output initiator file: ' 
read(*,'(A)') mame4 
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open 

(UNIT=14,FILE=fiiame4,STATUS='NEW',IOSTAT=FLERR4 
ERR=104) 

write(*,*) 'Writing data. . . ' 

END IF 
x=0 
i=0 
3=0 

do 22 i=l,LAMDA 
x= (i-l)+200 

call locate (xx2,NPTS,x j) 
if (0-eq.0).OR. O'.eq.NPTS)) then 
INTERV2=0 

else 

INTERV2= ((yy20+l)-yy20)) / (xx2(j+l) - xx2 
0))) * (x-xx2 0)) 
+ +yy2 (j) 

end if 

yi2 (i)=INTERV2 

IF ((INTMED.EQ.'Y') .OR. (INTMED.EQ. 'y ')) THEN 

if (O'.eq.O) .OR. O'.eq.NPTS)) then 
GO TO 22 

else 

write(14,'(lx,f7.2,i6,3fl2.2)') x, j, xx2 

0),xx2 0+1), 

+ INTERV2 
endif 
END IF 
22 continue 

IF ((INTMED.EQ. 'Y').OR (INTMED.EQ. 'y')) THEN 

close (UNIT=14) 

ENDIF 
ELSE 

ITER=0 
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do 24 ITER= 1 ,L AMD A 

yi2 (ITER)=1 
24 continue 



ENDIF 

write(*,*) 'Do you wish to process a light source file? (Y/N)' 
read(*,'(A)') LIGHT 

IF ((LIGHT.EQ.'Y').OR. (LIGHT.EQ.'y')) THEN 
ITER=0 

do 25ITER=1,NMAX 
xx3 (ITER) =0 
yy3 (ITER) =0 
yi3 (ITER)=0 
25 continue 

write(*,*) 'Enter the name of the light source file:' 
read(*,' (A)') fhame5 

open 

(UNIT=15,FILE=fhame5,STATUS='OLD',IOSTAT=FLERR5, 
ERR=105) 

ROWS=0 

do while ((EOF3.EQ.0) .AND. (ROWS.LT.NMAX)) 
ROWS=ROWS+l 

read(15,*,IOSTAT=EOF3) xx3 (ROWS), yy3 (ROWS) 

end do 

close (UNIT=15) 
NPTS=0 

IF (EOF3.NE.0) THEN 
NPTS=ROWS-l 

write(*,l (14,1X,A12)1)NPTS, 'points read.' 

ELSE 

NPTS=ROWS 

write (*,' (A28,I4,A12)') 'Too many data points! First ', 
+ NMAX,' points read. . . ' 

ENDIF 

call norm (yy3,NPTS) 
IF ((INTMED.EQ.'Y').OR. (INTMED.EQ. 'y ')) THEN 

write(*,*) 'Enter the name of the light source output file:' 
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read(V (A)') fhame6 
open 

(UNIT= 1 6,FILE=fiiame6,STATUS='NEW',IOSTAT=FLERR6 
ERR- 1 06) 

write(* *) 'Writing data. . . ' 

END IF 

x=0 
i=0 
j=0 

do 32 i=l,LAMDA 
x= (i-l)+200 

call locate (xx3,NPTS,x,j) 
if (O'.eq.O) .OR. (j.eq.NPTS)) then 
INTERV3=0 

else 

INTERV3= ((yy3(j+l) - yy3(j)) / (xx3(j+l) - 
xx3(j))) * (x-xx3 0)) 

end if 

yi3 (i) =INTERV3 

IF ((INTMED.EQ. 'Y') .OR. (INTMED.EQ. 'y ')) THEN 

if (G'.eq.O).OR. O'.eq.NPTS)) then 
GO TO 32 

else 

write(16,'(lx,f7.2,i6,3fl2.2)') xj,xx3 (j),xx3 (j-1), 

INTERV3 

endif 
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ENDIF 

32 continue 

IF ((INTMED.EQ. 'Y') .OR. (INTMED.EQ. 'y')) THEN 

close (UNIT=16) 

ENDIF 

ELSE 
ITER=0 

do 34 ITER= 1 ,L AMD A 

yi3 (ITER) -1 
34 continue 



30 



END IF 
ITER=0 

DO 40 ITER=1,LAMDA 

yc (ITER)-O 
40 CONTnSUE 

DO 55 i=l,LAMDA 

yc (i)-yil (i)*yi2 (i)*yi3 (i) 
55 CONTINUE 



write(*,*) 'Enter the filename for cumulative data-' 
read (*,'(A)') fhame7 
open 

(UNIT=17,FILE=mame7,STATUS-'NEW , ,IOSTAT=FLERR7 
ERR=107) 

write(*,*) 'Writing data. . . ' 
CALL integ (yc,LAMDA,area) 

write(*,'(lX,A26,All,F12.6)') "The area under the product', 
+'curve is:',area 

write(17,'(lX,A26,Al 1JF12.6)') 'The area under the product', 
+' curve is: ',area 

DO 60 i=l,LAMDA 
x= (i-l)+200 

write(17,' (1X,F6.1,2F11.2,F11.4,E15.6)') x,yil (i),yi2 

(i), 

+ yi3 (i) ,yc (i) 

60 CONTINUE 

close (UNIT= 17) 

101 IF (FLERR1 .NE. 0) THEN 
write(*,*) 'Unable to open substrate file!' 

END IF 

1 02 IF (FLERR2 .NE. 0) THEN 

write(*,*) 'Unable to create substrate output file!' 
END IF 

1 03 IF (FLERR3 .NE. 0) THEN 
write(*,*) 'Unable to open initiator file! ' 
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END IF 



1 04 IF (FLERR4 .NE. 0) THEN 

write(*,*) 'Unable to create initiator output file!' 
END IF 

1 05 IF (FLERR6 .NE. 0) THEN 

write(*,*) 'Unable to open light source file!' 
END IF 

106 IF (FLERR6 .NE. 0) THEN 

write(*,*) 'Unable to create light source output file!' 
END IF 

107 IF (FLERR7 .NE. 0) THEN 

write(*,*) 'Unable to create cumulative output file! ' 
END IF 

write(*,*) 'Program exiting normally. . . ' 
END 

SUBROUTINE locate (xx,n,x,j) 

INTEGER j,n 

REAL x,xx (n) 

INTEGER jljmju 

J1=0 

ju=n+l 

10 if(ju-jl-gt.l)then 
jm=0'u+jl)/2 

if ((xx (n).ge.xx (1)) .eqv. (x.ge.xx 0m))) then 
jl=jm 

else 

ju=jm 

endif 
goto 10 
endif 

if (x.eq.xx (l))then 
j=l 

else if (x.eq.xx(n))then 
j=n-l 

else 
endif 

return END 

SUBROUTINE baseln (yy,N) 
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INTEGER N,i 

REAL yy (N), minno, temp 

minno=yy (1) 
i=0 

temp=0 

DO 10i=2,N 

IF (yy (i) XT. minno) THEN 
minno=yy (i) 

END IF 
10 CONTINUE 

i=0 

DO 20i=l,N 

temp=yy (i)-minno 

yy (i)=temp 
20 CONTINUE 

END 

SUBROUTENffi norm (yy,N) 

INTEGER N,i 

REAL yy (N),maxno,temp 

maxno=yy (1) 
i=0 

temp=0 

DO 5 i=2,N 

IF (yy(i).GT.maxno) THEN 
maxno=yy (i) 

END IF 
5 CONTINUE 

i=0 

DO 10, i=l,N 

temp=yy (i)/maxno 

yy (i)=temp 
10 CONTINUE 

END 
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SUBROUTINE integ (yy,N,area) 

INTEGER N,i 

REAL yy (N),sum,area 

i=0 

sum=0 
area=0 

DO 10i=l,N-l 

sum=sum+ (yy(i)+yy (i+l))*0.5 
10 CONTINUE 

area=sum/ 100000 

END 



34 



APPENDIX 2 



Program to create uniformly spaced csv data from unevenly spaced 
tabular data 

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 

#define NMAX 3501 
#defme STRMAX 151 
#define FNMAX 81 
#define OUTPTS 801 

void locate(float xx[], unsigned long n, float x, unsigned long *j); 
void norm(float xx[] ? unsigned long int n); 
void baseline(float xx[], unsigned long int n); 

int mainO 

{ 

char fiiamein(FNMAX], string[STRMAX], *str, ptr, 
fiiameout[FNMAX], another; 

float xdata(NMAX], ydata[NMAX], xinter[OUTPTS+l], 

yinter[OUTPTS+l]; 

unsigned long int index, i, j; 

int choice; 

FILE *fpin, *fpout; 

another = 'Y'; 

do { 

for (i = 0; i <= NMAX-1 ; i++) { 
xdata [i] = 0; 

ydata [i] = 0; 

}; 

for (i = 0; i <= OUTPTS; i++) { 
xinter [i] = 0; 

yinter (i] = 0; 

}; 

printf("Enter name of the input file (80 chars max, no spaces): "); 
scanf("%s", fhamein); 

printf("File name is %s\n",fhamein); 

fpin = fopen(fnamein,"r"); 

if(fpin = NULL) { 
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printf("Cannot open %s\n",fhamein); 
exit(l); 

}; 

index = 1 ; 

while (1) { 

str_ptr = fgets(string,STRMAX-l,fpin); 
if(str_ptr = NULL) 

break; 
if (index = NMAX) 

break; 

sscanf (string, "%f %f , &xdata [index] , &ydata [index]); 
index++; 

}; 

f close (fpin); 

if((index = NMAX) && (strjrtr != NULL)) { 
index-; 

printf("Too many data points! Using first %d points 
only...\n",index); 
} 

else { 

index-; 

printf( M %d points read... \n", index); 

} ' 

printf("\nEnter option for data processing\n"); 
printf("l : Normalize the data after interpolation^"); 
printf("2: Baseline the data after interpolation^"); 
printf("3: First interpolate, then baseline and finally "); 
printf("normalize the data\n"); 
printf("4: Simply interpolate the data\n"); 
printf("5: Simply normalize the data\n"); 
printf("6: Simply baseline the data\n"); 
printf("or\n"); 

printf("0: to exit the program without any data processing\n"); 

printfC'XnEnter option (0-6): "); 

scanf("%d",&choice); 

if (choice = 0) 
exit(2); 

printf("\nEnter name of the output file (80 chars max, no spaces)' "V 
scanf("%s",fhameout); 

printf("File name is %s\n",fhameout); 
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fpout = fopen(rhameout,"w"); 

if (fpout = NULL) { 

printf("Cannot open %s\n",rhameout); 
exit(3); 

}; 

for(i = 1; i <= OUTPTS; i++) { 
xinterfi] = 200+((float)i-l); 
locate(xdata,index,xinter[i],&j); 
if «J = 0)110'= index)) ' 
yinter[i] = 0; 

else 

yinter [i] = (xinter [i] -xdata 0 ]) * ((ydata (j+1] -ydata [j ]) / 
(xdata [j+1] -xdata [j ])) +ydata [j ] ; 

}; 



if ((choice = 2) \\ (choice = 3)) { 
baseline(yinter,OUTPTS); 

}; 

if (choice = 6) { 

baseline(ydata,index); 

}; 

if ((choice = 1) || (choice = 3)) { 
norm(yinter,OUTPTS); 

}; 

if (choice = 5) { 

norm(ydata,index); 

}; 

if ((choice >= 1) && (choice <= 4)) { 
for (i = 1 ; i <= OUTPTS-1 ; i++) { 
fprintf(fpout,"%13.5E, ",yinter[i]); 

Q)rintf(fpout,"%13.5E\n",yinter[OUTPTS]); 
else 

if ((choice == 5) || (choice = 6)) { 
for (i = 1 ; i <= index- 1; i++) 
fprintf(rpout,"% 1 3 ,5E, n ,ydata[i]); 

} 5 

rprintf(fpout,"%13.5E\n",ydata[index]); 

} y 

fclose(fpout); 
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printff'File %s written.\n\n",fhameout); 
prints-Process another file (Y/y/N/n)?: "); 
scanf("%ls",&another); 
} while (another = T \\ another = 'y'); 

printf("Exiting...\n"); 
retum(O); 

} 

void locate(float xx[], unsigned long n, float x, unsigned long *j) 

{ 

unsigned long jujmjl; 
int ascnd; 

ji=0; 

ju=n-H; 

ascnd==(xx[n] >= xx[l]); 
while (ju-jl > 1) { 

jm=(ju+j 1)» 1; 

if (x >= xx [j m] = ascnd) 

else 

ju=jm; 

} 

if(x = xx[l]) 
else if (x = xx[n]) 

*j=n-i; 

else 

*j=ji; 

} 

void norm(float xx[], unsigned long int n) 

{ 

unsigned long int i; 
float maxdata, temp; 

maxdata = xx [1] ; 
temp = 0; 

for(i = 2;i<=n; i++) { 
if(xx[i] > maxdata) 
maxdata = xx[i]; 

}; 

for(i = l;i<=n; i++) { 
temp = xx[i]/maxdata; 
xx [i] = temp; 

}; 
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} 



void baseline(float xx[], unsigned long int n) 

unsigned long int i; 
float mindata, temp; 
mindata = xx [ 1 ] ; 
temp = 0; 

for(i = 2; i <= n; i++) { 
if(xx[i] < mindata) 
mindata = xx [i] ; 

}; 



15 for(i = l;i<=n;i++) { 

temp = xx(i] - mindata; 
xx(i] = temp; 

}; 

} 
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APPENDIX 3 

Program for Determining Strength of Wavelength Response in a 

Region 

#include 

<stdio.h> 

#include 

<stdlib.h> 

#include 

<math.h> 

#define NMAX 
3501 

#defme STRMAX 
151 

#define FN MAX 
81 

#define OUTPTS 
801 

void locate(float xxQ, unsigned long n, float x, 
unsigned long *j); 

void norm(float xxQ, unsigned long 
int n); 

void baseline(float xxQ, unsigned 
long int n); 

void partinteg (float xxQ, unsigned long int x1 , unsigned long int 

x2, 

float *area); 
int 

main() 
{ 

char fnamein[FNMAX], string[STRMAX], *str_ptr, 
fnameout[FNMAX], another; 

float xdata[NMAX], ydatafNMAX], xinter[OUTPTS+1], 
yinter[OUTPTS+1]; 

float 

totalarea.aA.aB.aC.aD.aE.aF.aG.aH.al.aJ.a 
K; 

unsigned long int index, i, 

j; 

int choice; 
FILE *fpin, 
*fpout; 
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another = 'Y'; 

printf("Contact Rajdeep S. Kalgutkar, SRC-CRC 7-3003, for 
further info\n"); 

do{ 

for (i = 0;i<= NMAX-1; 

xdata[i]=0; 
ydata[i]=0; 

}; 

for (i = 0; i <= OUTPTS; 
i++){ 

xinter[i]=0; 
yinter[i]=0; 

}; 

printf("\n Enter name of the input file (80 chars max, no 
spaces): "); 

scanf("%s",fnamein); 

printf("File name is 
%s\n",fnamein); 

fpin = 
fopen(fnamein,"r"); 

if (fpin -- 
NULL) { 

printffCannot open %s. 
Exiting... \n",fnamein); 

exit(l); 

}; 

index = 1 ; 

while (1){ 
str_ptr = fgets(string,STRMAX- 
1 ,fpin); 

if(str_ptr == NULL) 

break; 
if(index == NMAX) 

break; 
sscanf(string,"%f 
%f,&xdata[index],&ydata[index]); 
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index++; 

}; 

fclose(fpin); 

if((index == NMAX) && (str_ptr != NULL)) { 
index—; 

printffToo many data points! Using first %d points 
on!y...\n , \index); 
} 

else { 
index—; 

printf("%d points 
read...\n",index); 

}; 

printfC'VnEnter option for data 
H; processing\n M ); 
S printf("1 : Simply interpolate the 

2 data\n"); 

printf( M 2: Normalize the data after 
y interpolation^"); 
yp printf("3: Baseline the data after 

O interpolation^"); 

f printf("4: First interpolate, then baseline and finally 

N- »\, 

51 printf("normalize the 

m data\n"); 

O printf("or\n"); 

H printf("0: to exit the program without any data 

processing\n"); 

printffViEnter option (0- 
4):"); 

scanf("%d M ,&choice); 

if (choice == 0) 

exit(2); 

printf("\nEnter name of the output file (80 chars max, no 
spaces): "); 

scanf( H %s",fnameout); 

printf("File name is 
%s\n",fnameout); 

fpout = 

42 



fopen(fnameout,"w"); 

if (fpout == 
NULL){ 

printffCannot open %s. 
Exiting. ..\n",fnameout); 

exit(3); 

}; 

for(i = 1; i <= OUTPTS; 
i ++ ) { 

xinterfi] = 200+((float)i- 

1); 

locate(xdata,index,xinter[i],&j); 

if((j == 0) || G == 
index)) 

yinter[i] = 0; 

else 

yinter[i]=(xinter[i]-xdataD])*((ydataD+1]- 
ydataD])/ " ' 

(xdata[j+1]-xdataD]))+ydata[j]; 

if ((choice == 3) || (choice == 4)) { 
base!ine(yinter,OUTPTS); 

}; 

if ((choice == 2) || (choice == 4)) { 
norm(yinter,OUTPTS); 

}; 



partinteg(yinter,51 ) OUTPTS,&totala 
rea); 



partinteg(yinter,51,101,&a 
A); 

partinteg(yinter, 1 01 , 1 51 ,& 
aB); 

partinteg(yinter, 1 51 ,201 ,& 
aC); 
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partinteg(yinter,201 ,251 ,& 
aD); 

partinteg(yinter,251 ,301 ,& 
aE); 

partinteg(yinter,301 ,351 ,& 
aF); 

partinteg(yinter,351 ,401 ,& 
aG); 

partinteg(yinter,401 ,451 ,& 
aH); 

partinteg(yinter,451 ,501 ,& 
al); 

partinteg(yinter,501 ,551 ,& 
aJ); 

partinteg(yinter,551,OUTPTS,&aK); 

fprintf(fpout,"The total area is: 
%14.6E\n",totalarea); 

fprintf(fpout,"The area under region A is: 
%6.2f%%\n",aA*1 00/totalarea); 

fprintf(fpout,"The area under region B is: 
%6.2f%%\n",aB*1 00/totalarea); 

fprintf(fpout,"The area under region C is: 
%6.2f%%\n",aC*1 00/totalarea); 

fprintf(fpout,"The area under region D is: 
%6.2f%%\n",aD*1 00/totalarea); 

fprintf(fpout,"The area under region E is: 
%6.2f%%\n",aE*1 00/totalarea); 

fprintf(fpout,"The area under region F is: 
%6.2f%%\n",aF*1 00/totalarea); 

fprintf(fpout,"The area under region G is: 
%6.2f%%\n",aG*1 00/totalarea); 

fprintf(fpout,"The area under region H is: 
%6.2f%%\n",aH*1 00/totalarea); 

fprintf(fpout,"The area under region I is: 
%6.2f%%\n",al*1 00/totalarea); 

fprintf(fpout,"The area under region J is: 
%6.2f%%\n",aJ*1 00/totalarea); 

fprintf(fpout,"The area under region K is: 
%6.2f%%\n\n",aK*1 00/totalarea); 
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for (i = 1;i<= OUTPTS-1; { 
fprintf(fpout,"%1 3.5E, ",yinter[i]); 

}; 

fprintf(fpout,"%13.5E\n",yinter[i]); 
fclose(fpout); 

printf("File %s 
written.\n\n",fnameout); 

printf("Process another file 
(Y/y/N/n)?: "); 

scanf("%1 s",&another); 
} while (another == *Y' || another 
== V); 



printf("Exiting...\n 

"); 

return(0 

); 

} 

void locate(float xxQ, unsigned long n, float x, 
unsigned long *j) 

{ 

unsigned long jujmjl; 
int ascnd; 

jl=0; 

ju=n+1; 
ascnd=(xx[n] >= xx[1]); 
while Ou-jl > 1) { 
jm=0u+jl) » 1; 
if (x >= xx[jm] == ascnd) 

jl=jm; 
else 

ju=jm; 
} 

if (x ==xx[1]) 

1=1; 
else if(x == 
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xx[n]) 
*j=n- 

1; 

else 

*j=ji; 

} 

void norm(float xxg, unsigned long 
int n) 

{ 

unsigned long 
int i; 

float maxdata, 
temp; 

maxdata = 
xx[1]; 
temp = 

0; 

for(i = 2; i <= n; i++) { 
if(xx[i] > 
maxdata) 

maxdata = 
xx[i]; 

}; 

for(i = 1 ; i <= n; i++) { 
temp = xx[i]/maxdata; 
xx[i] = temp; 

}; 

} 

void baseline(float xxQ, unsigned 
long int n) 

{ 

unsigned long 
int i; 

float mindata, 
temp; 

mindata = xx[1]; 
temp = 

0; 

for(i = 2; i <= n; i++) { 
if(xx[i] < 



mindata) 

mindata = 
xx[i]; 

}; 

for(i = 1; i <= n; i++) { 
temp = xx[i] - mindata; 
xx[i] = temp; 

}; 

} 

void partinteg(float xxQ, unsigned long int x1, unsigned long int 
x2, 
float *area) 

{ 

unsigned long 
int i; 

float temp; 
temp = 



for(i = x1 ; i <= x2 - 1 ; i++) 
temp = temp + (xx[i] + xx[i+1])/2; 

}; 

*area = temp; 
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APPENDIX 4 



SRC Curing Resource dB 4 Query Select2 

5 Sub Initialize 

Dim ses ses As New NotesSession 
Dim db_db As NotesDatabase 
Dim view view As Notes View 

Dim note jaotel As NotesDocument, note_note2 As NotesDocument 
10 Dim i cnt As Integer, i add As Integer 

Set db db =ses ses.CurrentDatabase 

Set note notel=ses ses.DocumentContext 

1 5 Redim Preserve arr WavelengthRegion(O) air WavelengthRegion(O) _ " c 

If note_notel.Selectionl(0) o Or rwte_noteLSelection2(0) o Then 
U If note noteLSelectionl(O) o Then 

O Set view_yiew = db_db.GetView("By NotelD") 

p20 

> If note notel.Selection2(0) o "'Then 

Set notejiote2^iewj^iew.GetDocumentByKey(Right("00000000" & 
2 note_notel.Selection2(01, 8)) 

Q E^ se 

a 25 Set note note2 = view_view.GetDocumentByKey(Right("00000000" & 
H notejQOtel.Selectianl(0), $)) 

O End If 

:4 If Not (notejiote21s Nothing) Then 

tj30 If note note2.Hasltem('WavelengthRegion") Then 
i_cnt = -1 

Forall vals In note_note2.WavelengthRegion 
Ifvalso""Then 

i_cnt = i_cnt + 1 

3 5 Redim Preserve arr WavelengthRegion(i cnt) 

arr_WavelengthRegion(i cnt) = vals 

End If 
End Forall 

40 End If 

End If End If 

If note notel.Type(O) = "S" Then 

Set view- View = db db.GetView("Substrate") 
45 Elseif note _notel.Type(0) = "P" Then 

Set view- view = db_db. Get Vie w('Tnitiator Sensitizer") Else 
Set view- view = db_db.GetView("LightSource") End If 
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' Set note note2 =view-view. GetFirstDocument 
i_cnt = -1 

Do While Not (note note2 Is Nothing) 
If note_note2.Name(0) o Then 
L add = True 

If air WavelengthRegion(O) o Then 
Ladd = False 

Forall valsl In note _note2.WavelengthRegion 
Forall vals2 In arrJWavelengthRegion 
If valsl = vals2 Then 
i__add = True 
Exit Forall 

End If 
End Forall 

If Ladd Then 

Exit Forall 

End If 
End Forall 
End If 

If Ladd Then 

i cnt = i cnt + 1 

Redim Preserve arr_names(i cnt) 

arr names(i_cnt) = note_note2.Name(0) 

End If 

End It 

Set note note2 = view- view. GetNextDocument(note_note2) 

Loop 

note __notel .Names = arr_names 

End Sub 

SRC Curing Resource dB 4 Query Select2 Save Agent 

Sub Initialize 

Dim ses sesAs New NotesSession 
Dim db_db As NotesDatabase 
Dim view view As Notes View 

Dim note jiotel As NotesDocument, note_note2As NotesDocument 

Set db_db = ses ses.CurrentDatabase 
Set note note) = ses ses.DocumentContext 
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Select Case notenotel.Type(O) 
Case "S" 

Set view view = db_db.GetView("(Substrate)") 
5 Set note_note2=view view.GetDocumentByKey(note_notel.Substrate(0)) 
Case "P" 

Set view-view db-db.GetView("(InitiatorSensitizer)") 
10 Set note note2 -view view.GetDocumentByKey(note notel.Photolnitiator(O)) 
Case "L" 

Set view-view =db db.GetView("(LightSource)") 
15 Set note_note2 = view_view.GetDocumentByKey(note_notel.LightSource(0)) 
End Select 

g If note_notel.MexWction(0) = "Add" Then 

[if 20 If note_notel.Selectionl (0) <> "" Then 

*?■ Print "[!" + note notel.dbname(O) + "/QuerySelectionl?OpenForm&" & 

hi note_notel.Selectionl (0) & "&" & note_note2.Noteld & ")" 

§ Else 

O Print "[!" + note notel.dbname(O) + "/QuerySelectionl?OpenForm&" & 

f 25 note_note2.Noteld & "]" 

JT ; End If Elseif note notel.NextAction(O) = "Separate" Then 

% If note_notel.Selectionl(0) o "" Then 

pi Print " [/" + note notel.dbname(O) + "/QuerySelectionl?OpenForm&" & 

11; note_notel.Selectionl(0) & "&" & note note2.Noteld & ")" 

S30 Else 

Print "[!" + note notel.dbname(O) + "/QuerySelectionl?OpenForm&" & 
note_note2.Noteld & ")" 
End If 

Elseif note _notel.NextAction(0) = "Separate" Then 
35 If note _notel.Selection2(0) o "" Then 

Print "[/" + note notel.dbname(O) + "/QuerySelectionResults?OpenForm&" & 

note_notel. Selection) (0) & "&" & note note).Selection2(0) & _ 

"&" & note_note2.Noteld & "]" 

Elseif note notel.Selectionl(O) o "" Then 
40 Print- "[" + note notel.dbname(O) + "/QuerySelectionResults?OpenForm&" & 

note_notel.Selectionl (0) & "&" & note note2.Noteld & ")" 

Else 

Print "[/" + note_notel.dbname(0) + "lQuerySelectionResults?OpenForm&" & 
note_noteMoteld & "]" 
45 End If 
Else 

If note notel.Selection2(0) o "" Then 
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Print "[/" + note_notel.dbname(0) + "/QuerySelectionOverlayResults?OpenForm&" & 

notenotel.Selectionl(O) & "&" & note_notel.Selection2(0) & 

"&" & note_note2.Noteld & ")" 

Elseif note_notel.Selectionl(0) o ""Then 

Print "[/" + note_notel.dbname(0) + "/QuerySelectionOverlayResults?OpenForm&" 

& note_notel. Selection) (0) & "&" & note note2.Moteld & ")" 

Else 

Print "[/" + note_notel.dbname(0) + "lQuerySelectionOverlayResults?OpenForm&" 
& note_note2.Noteld & ")" 
End If 

End If 
End Sub 
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SRC Curing Resource dB 4 Query Overlay Open Agent 

Sub Initialize 

pim ses_ses As New NotesSession 
Dim db db As NotesDatabase 
Dim view- view As NotesView 

Dim note notel As NotesDocument. note _note2As NotesDocument 
Dim icntAs Integer, i addAs Integer 

Set db db = ses ses.CurrentDatabase 
Set note notel =ses ses.DocumentContext 

S et view- view = db_db. G etView("By N otel D ") 

If note notel. Selectionl(O) o Then 

S et note-note2 = view- view. G etD ocumentByKey(R fight("00000000" + note_noteL 
S election"! (0). 8)) 

If Not (note note2 Is Nothing) Then 

note _notel.datal = note_note2.EmissData 

note notel.maxfreql = note_note2.MaxFreq 
End If End If 

If note_notel.Selection2(0) o Then 

Set note-note2 = view view.GetDocumentByKey(Right("00000000" + 
note_notel.Selection2(0), 8)) 

If Not (note note2 Is Nothing) Then 

note notel.data2=note_jiote2.EmissData 

note _notel.maxfreq2 = note_note2.MaxFreq 
End If End If 

If note notel.Selection3(0) o Then 

Set note-note2 = view_view.GetDocumentByKey(RightC t 00000000 ,, + 
notejnotel.Selection3(0), 8)) 

If Not (note _note21s Nothing) Then 

note notel.data3 = note note2.EmissData 
note jnotel.maxfreq3=note note2.MaxFreq 
EndlfEndlfEndSub 
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import java,awt.*; import java.awt event.*; import java.applet*; 

public class SRC Charts extends Applet { int gi count; 

double GetHMax(String strjn) { String str_current; double dbl hmax; 

str_current = ""; dblhmax = 0; for Tint i cnt = 0; i cnt < str_in.length(); i cnt++) { 

if(stx in. region Match es(i cnt, 0, 2)) 

if(Double.valueOf(str current),doubleValue() > dblhmax) 

dbl hmax = Double. valueOf(str current).doubleValue(); 

str current = ; 

i_cnt-H-; 

gi count-H-; ) else { 

str current = str current.concat(str in.substring(i cnt, i cnt + 1)); 
) 

if(str_current.length() > 0) { 

if(Double.valueOf(str_current).doubleValueO > dbl_hmax) 
dbl hmax = Double.valueOf(str_current).doubleValue(); gi count-H-; ) return dbl hmax; 



int StringTolnt(String strain, double dbl hmax) { double dbl_pos; 

dbl_pos = getSize().height - (25 + (Double.valueOf(str_in).doubleValueQ * 
((getSizeO-height - 50) / dbl hmax))); return (int)dbl-pos; ) 

void drawChartLine(Graphics g, String str_in, String strjype, double dbljmaxfreq) { 
double dbl_x, dbl_inc, dbljhmax; String str last, str next; 

str last strjnext = ; dbl x = 25; gi count = 0; 

dbl_hmax = GetHMax(str_in); if(str type, equals IgnoreCase("S")) dbl hmax = 100; 
dbl inc = ((((double)getSizeQ.width - 50) / gi count) * ((dbl maxfreq - 200) / 800)); 
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for (int i cnt = 0; i_cnt < str_in.length(); ) i c { 

nt++)if(str_in.regionMatches(i cnt, 0, 
2))if(str_last.length(> 0) { 

{ 

g.drawLine((int)dbl x, StringTolnt(str_last, dbl hmax), (int)(dbl x + dbl inc), 
StringTolnt(str_next, dbl hmax)); 
dbl x = dbl x + dbl inc; 

str last = str_next; str next = ""; i cnt++; 

else { str next = str next.concat(str in.substring(i cnt, i cnt + 1)); 
) 

if(str_next.length() > 0) 

g.drawLine((int)dbl x, StringTolnt(str_last, dbl hmax), (int)(dbl x + dbl-inc) 
StringTolnt(stmext, dbl_hmax)); ) 

public void paint(Graphics g) { double dbl x, dbl-y; 

g.setColor(Color.black); g.drawLine(0, 0, getSize().width, 0); g.drawLine(25, 
getSize().height - 25, getSize().width - 25, getSize(). height - 25); g.drawLine(25, 25 
25, getSize().height - 25); 

for(int i cnt = 0; i_cnt < 9; i cnt++) { 

dbl_x = 25 + ((double)i cnt * (((double)getSize().width - 50) / 8)); 
g.drawLine((int)dbl x, getSize(). height - 25, (int)dbl x, getSize().height - 20); 
g.drawString(String.valueOf((i cnt * 100) + 200), (int)dbl x - 8, getSize(). height - 5); ) 

for(int i cnt = 0; i_cnt < 1 1; i_cnt++) { 

dbl-y = 25 + ((double)i cnt * (((double)getSize().height - 50) / 10)); 
g.drawLine(20, (int)dbl-y, 25, (int)dbl-y); 

g.drawString(String.value0f(100 - (i cnt * 10)), 1, (int)dbl-y + 5); ) 
g.drawString("Data Overlay", (getSize().width / 2) - 30, 12); 

g.setColor(Color.red); drawChartLine(g, getParameter("Data 1 "), getParameter("Type 
1 "), Double.valueOf(getParameter("Max Freq 1 ")). double Value()); 
g.setColor(Color.blue); drawChartLine(g, getParameter("Data 2"), ' 
getParameter("Type 2"), Double. valueOf(getParameter("Max Freq 
2")).doubleValue()); g.setColor(Color.green); drawChartLine(g, getParameter("Data 
3"), getParameter("Type 3"), Double.valueOf(getParameter("Max Freq 
3")).doubleValue()); ) ) 
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